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1. Introduction 

The dimensional regularization with the minimal subtraction has become everyone's 
favorite method for perturbative calculations of renormalizable field theories. [|l|] The 
method has three advantages. First, calculations are simple. Secondly, it is mass indepen- 
dent. Thirdly, the singularities higher than 1/e are determined completely by the simple 
1/e poles. § 

Consider the ((/>^)4 theory as an example. The theory 

has three parameters g^, m^, and A with scaling dimensions 4, 2, and 0, respectively. 
These parameters are conjugate to the operators 1, ^^/2, and (p^/Al whose scaling dimen- 
sions are 0, 2, and 4, respectively. Under the renormalization group (RG), the coordinate 
distance r transforms into re~', while the renormalization point is always fixed at r = 1. 
Hence, the RG acts toward the infrared limit, and it differs from the standard definition 
by rescaling. Then, under the RG the parameters transform as follows: 

I = {2 + PU^))m' (1-1) 
al 

^A = /3.(A). 

As we can see, the scaling dimensions of the parameters are additively preserved under 
the RG; only terms of scaling dimension 4 are allowed in the first equation, while only 
terms of scaling dimension 2, are allowed in the second and third, respectively. In fact 
in ref. it was shown that we can always choose parameters such that the scaling di- 
mensions of the parameters are additively preserved in any field theory with an ultraviolet 
fixed point. These parameters were also shown to be appropriate to describe the short 
distance physics. Locality of the theory implies that the RG equations involve only inte- 
gral powers of these parameters. Hence, the RG equations become finite polynomials of 
dimensionful parameters whose coefficients are power series in dimensionless parameters. 
Within perturbation theory the results of ref. [§] applies also to the (p'^ theory. The min- 
imal subtraction with dimensional regularization provides an example of such a choice of 
parameters. The structure of ( PTT] ) is called mass independent, since the beta functions 
depend only on the dimensionless parameter A. 

Now, in the minimal subtraction scheme in the dimensional regularization, the beta 
functions and the anomalous dimensions of composite operators are directly related to the 
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simple 1/e poles of the unrenormalized correlation functions. In other words, nontrivial 
renormalization properties of the theory, i.e., nonvanishing beta functions and anomalous 
dimensions, demand that the bare correlation functions have simple 1/e poles. Since the 
higher order poles in 1/e are related to the simple 1/e poles, we can say that all divergences 
in e are expected consequences of nonvanishing beta functions and anomalous dimensions. 

So far we have explained the advantages of the dimensional regularization with the 
minimal subtraction. We should not be totally happy with this method, however, since 
the physical meaning of the divergences in e is unclear. The purpose of this paper is to 
consider physical singularities of correlation functions at short distances and relate the 
short distance singularities, rather than the unphysical singularities in e, to the anomalous 
dimensions. Our results are valid beyond perturbation theory, and we will use QCD in 
four dimensional euclidean space as an example. 

The paper is organized as follows. We summarize the relevant facts on the RG equa- 
tions in sect. 2 and those on the operator product expansions in sect. 3. Then, we introduce 
the main formula in this paper in sect. 4 that describes the change of correlation functions 
of composite operators under an infinitesimal change of parameters. In sect. 5 we derive a 
relation between the operator product coefficients and anomalous dimensions of composite 
operators by considering the consistency between the main formula of sect. 4 and the RG 
eqs. In sect. 6 we make yet another consistency check of the main formula of sect. 4, com- 
ing from commutativity of derivatives. In sect. 7 we introduce a convention for composite 
operators. We give concluding remarks in sect. 8. 

2. Renormalization group equations 

In this section we introduce relevant facts on the RG equations in QCD with massive 
quarks in four dimensional euclidean space. The theory is characterized by three parame- 
ters g-L, m, and qe- The parameter (/i is an additive constant to the lagrangian density, and 
its scaling dimension is four. The parameter m is the quark mass parameter with scaling 
dimension one, and is the strong fine structure constant with scaling dimension zero.0 
We denote the operators conjugate to these parameters by 1, Om, and Oe, respectively. 
The operator 1 is the identity operator. The operator Om is the mass density operator 

^ We do not consider a gauge fixing parameter, since we will only consider gauge invariant 
operators in this paper. 
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'i/j^fj, and Oe corresponds to the energy density operator F^^F^^. These three operators 
are the only gauge invariant scalar operators with scaling dimension less than or equal to 
four that conserve C and P. There is no other independent operator with these properties; 
d^ii^'j^i^) = and ip^^^D^ip = rm()i() by the eqs. of motion. 
We can write down the RG eqs. in the following form: 

Jl 91= ^91 + -^Pi{9e) 

^^m = il+Pmi9E))m (2.1) 
9e= l3Ei9E)- 

The above form is determined by the requirement that the scaling dimensions of the 
parameters are additively conserved under the RG. In ref. it was shown that we can 
always choose such parameters. In the first eq. in ( |2.1[ ) the coefficient of g-i is simply 4, 
since the parameter g-i, being an additive constant in the lagrangian, has nothing to do 
with interactions. Locality of the theory implies that the beta functions can be expanded 
in powers of qe near Qe = 0- 

l3i{gE) = /3i,o + l3i,i9E + ■■■ 

Pm{9E) = l3m,i9E + ^^9e + ■■■ (2.2) 
Pe{9e) = -^9e + -^9e + - ■ 
Let F{gx,'m, qe) be the free energy density. Then, by definition of the RG, we find 

= 4F. (2.3) 

Since 



we find, from (|2.3|) , and (|2.4|) , that the composite operators Om and Oe satisfy the 

RG eqs. 



/7 fyi 

-Om = (3 - f3m{9E))Om - -^Pi(9e)1 

-Oe = (4 - f3E'{9E))OE - m(3m'{9E)Om - -rrf3i'{9E)l 
dl 4! 



(2.5) 



In principle there can be additional scalar operators whose expectation values vanish, for 
example total derivative operators. But there is no other scalar operators of dimension 
less than or equal to four, and the above equations are exact. 
We now introduce RG eqs. for general composite operators. 

Let {Oi} be a basis of gauge invariant scalar operators, where Oi has scaling dimension 
Xi- Then the RG eqs. take the form 

j^O, = x,0^ + J]r,,,(m, (2.6) 

3 

where the anomalous dimensions j can be nonvanishing only if Xi — Xj is a non-negative 
integer. More precisely, the anomalous dimensions are finite polynomials of m with coef- 
ficients as power series in qe- 

^i,jirn,gE) = 7^,J(^E)m^'-^^ (2.7) 

where ■jij{gE) is a power series. We assume that 

r.,,(0,0) = (2.8) 

so that at the UV fixed point m = qe = there is no mixing of operators. In the matrix 
notation, in which Oi is an infinite dimensional column vector O, we can write eqs. ( |2.6| ) 

as 

j^O = XO + Tim,gE)0, (2.9) 

where Xa . diagonal matrix. 

Before closing this section, we will examine the RG eqs. (|2.1| ) further for later con- 
venience. We define the running parameters QEi^', 9E),'nT'{l;'m, qe), and gxil] g\,m, qe) as 
the solutions of the RG eqs. 

d _ 

9e{1;9e) = (^E{9E{k9E)) 

d 

— m{l] m, Qe) = (1 + l3m{9E{k 9E)))m{l] m, qe) (2.10) 

^ 9i\h9i:'m,9E) = Agx[l]gx,m,gE) H Pi{9e{1; 9e)), 

with the initial conditions 

9e{0]9e) = 9E: fh{0]m,gE) = m, gi{Q] gi,m, gE) = gi- (2.11) 
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We will suppress the initial conditions from now on, and we will denote, for example, 
simply gE{i) instead of gE{i',gE)- Then we note that 



5f£;(lnr), m(lnr), (7i(lnr) 



are all RG invariants, since the dependence on the change of the coordinate r under the 
RG is canceled by the dependence on the change of the initial parameters gi,m,gE under 
the RG. Using the beta functions we can write the running parameter rh{l) as follows: 



m{l) = me^E{gE{l),gE), 



(2.12) 



y 



where E{y,x) is defined by 

E{y,x) = exp 
Note that the function E satisfies 



d d 
(^E{y)-Q;^E{y, x) = (3m{y)E{y, x), (3E{x)—E{y, x) 



-E{y,x)(3^{x) 



(2.13) 



^2.14) 



3. Operator product expansions 

In this section we summarize relevant facts on the operator product expansions 
(OPE).[H We will be interested in the two types of operator products: 



O^(r)O(0) = Cm{r-.m,gE)0{Q) + o 



OE{r)0{Q) = CE{r- m, gE)0[^) + o 



(3.1) 



Here, both Cm and Ce are part of the operator coefficient functions that are at least as 
singular as l/r"*. In taking the operator products we have taken the average over the 
orientation of the coordinate vector r^; without the average the right-hand sides of ( |3.1J ) 
would include non-scalar composite operators. 

One of the most fundamental properties of OPE is that for a fixed coordinate r the 
coefficient functions can be expanded in powers of the small parameters m,gE- We will 



find that this analyticity, together with the RG, implies that the coefficient functions are 
finite polynomials in m whose coefficients are infinite power series in qe, i.e., 



n=0 



^ 777/ ^ 

n=0 



(3.2) 



We note that this implies that Cmij can be nonvanishing only for Xi > xj + 1, while CEij 
can be nonvanishing only for Xi > Xj. Let us derive the above results. 
The coefficient functions are closed under the RG, and we find 

^ Cm{r;m,gE) = - Pm{gE))Cm{r;m, qe) + [X + T{m, qe), Cm{r;m, qe)] 

4 CE{r; m, qe) = (4 - I3e {9e))C E{r; m, qe) - m(3rn' i9E)Cni{r; m, qe) ^^'^^ 
dl 

+ [X + r(m, gE),CE{r; m, qe)]- 
In order to solve these RG eqs. we introduce a matrix G{r;m, qe) that satisfies 

^^(r; m, qe) = {X + r(m, qe)) G{r; m, qe) (3.4) 
dl 

and the initial condition G{1; m, qe) = 1. The solution is given by 



G{r; m, Qe) = Texp 



dx 

gsilnr) (^e{x ^ 



T(m(lnr)i?(x, (/^(Inr)), x) 



(3.5) 



where T denotes the increasing ordering of x from right to left. Due to eq. ( |2.7| ) 
rij{m, qe) = rn^'~^^ JijigE), the matrix element Gij{r;m, qe) can be nonvanishing only 
if Xi — Xj is a non-negative integer. As far as the power of r is concerned (i.e., we ignore log- 
arithmic corrections), m(lnr) is proportional to r from eq. ( |2.12| ). Hence the dependence 
of Gij{r; m, qe) on the powers of m and r is given by 



Gij{r;m,gE) oc 



(3.6) 



Similarly, we find 



(3.7) 



It is also helpful to note that the matrix G{r] m, qe) satisfies 
d 



dhir 



G{r]m,gE) = -G{r]m, gE){X + F {m {In r) , gE{lnr))) . (3.^ 



Using the matrix G we can solve the RG eqs. ( |3.3| ) as 
follows: 

Cm{r; m, qe) = ^E{gE{lrir), gE)G{r; m, 5f£;)iJ^(m(ln r), (7^;(ln r))G'"\r; m, gE) 
CE{r]m,gE) = \ ^^^f^}^^^^^ G{r]m,gE)HE{rh{lnr),gE{lnr))G~'^{r]m,gE) (3.9) 

Tfl 

+ -5—^ — r(/3m((7E(lnr)) - P^{gE))Cm{r;m, gE), 
Pe{9e) 

where we define 

Hm{rn,gE) = Cm{l;m,gE), HE{m, gE) = CE{l;rn, gE). (3.10) 

Prom the analyticity assumption, both Hm{m, gE) and HEijn.gE) are power series in m 
and gE- Suppose H^ij{m, gE) has a term proportional to m"^. Then, from eq. ( |3.9|) 

we find that Cmi,j{f",'>^, Se) has a term proportional to rn^ /r^+^i~^i~'^. Since 
Cm{r; m, gE) must be at least as singular as by definition, we must have n < Xi—xj — l. 
Hence, Hmi,j{rn, gE) is a polynomial of degree Xi — Xj — 1 as far as m is concerned. The 
m"^ term in Hmij can also contribute to Cmi',j'{r;m, gE), where x[ > Xi > Xj > x'j. But 
from (|3.6| ) and ( |3.7| ) we find again that n < Xi ~ Xj — 1. Similarly, we can conclude that 
HEi,j{m, gE) is a polynomial of degree Xi — Xj with respect to m. To summarize, we find 
the structure 

H^,^j{m,gE)= Yl ^^^r^lj^SE), 

■ (3.11) 



HE^,J{m,gE)= V ^H^EiAgE)- 



n=0 



4. Variational formulas 

We will introduce a main formula that gives operator realization of the derivatives with 
respect to m and gE in this section. We will see that the consistency of this formula with 
the RG eqs. gives a relation between the operator coefficients and anomalous dimensions. 
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So far we have been saying casually that the parameter m and the operator Om are 
conjugate to each other. (The following discussion applies to qe and as well.) The 
precise meaning of the conjugacy relation is that the derivative — ^ is realized by the 



operator Om- One example is eq. (|2.4|) . This conjugacy relation must be valid not only for 



the simple expectation value as (|2.4| ) but also for any correlation functions of composite 



operators. So, we expect a formula of the type 
d 



gJO,Ari)...0,„{rn))m,g, = J {Om{r)0,Ari)...0,^{rn)rm,,,. (4.1) 
where 

{Omir)0,Ari)...OMrm,g, ^ {(Omir) - (O^)^,, JO., (ri)...0,„ (r^))^,,,. (4.2) 

We immediately notice, however, that the above formula is not well defined due to the short 
distance singularities in the product of Om and Oi. We emphasize that these short distance 
singularities are physical singularities that exist in QCD even after renormalization. (We 
will make a brief remark on the relation between these short distance singularities and the 
singularities we encounter in perturbative calculations in the concluding section.) There 
must be a way of fixing this problem in a local way, since the singularities result from the 
UV physics rather than the IR physics. Thus, we postulate the validity of the following 
variational formulas in QCD (an analogous formula was mentioned in ref. Q): 
d 



dm 



(0,,(ri)...0.„(rj)^,,, = {0*mO,,{n)...OM) 



m,gE 



m,gE 



^ hm [ / d\ {Om{r)0,Ari)...O^Sr^)) 

n 

mikji^i ^5 9e) {Qii {fl)---Oj{rk)...Oi^ {fn))m,gE 

k=l 

— (0,,(ri)...0,„(rj)^,,, = {0*EO,M)■■■O^Arn))^,gE 



(4.3) 



m,gE 



dgE 

^\im \ [ d'^r {OE{r)O.Ari)...OM) 

n 

+ AE^^,J{e■, m, Qe) {Oi^ {ri)...Oj{rk)...Oi^ (^n))m,gB 

k=l 



where we define 



Jl>r>£ 

AEi,j{e;m,gE) = - / d'^r C Ei,j{r;m, Qe) + CEi,j{m, Qe)- 

Jl>r>e 



(4.4) 
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The variational formulas (^4.3| ) explain why it is important to treat the operator coefficients 
at least as singular as l/r'^; less singular coefficients are irrelevant for the subtractions. The 
finite counterterms Cm, ce are necessary in order to compensate the arbitrariness involved 
in the subtracting procedure. The finite 

counterterms must have the same structure as the operator coefficients Cm, Ce (for 
a fixed r), and they can be expanded as 



T1 = 



For later convenience we define the maximal parts of the counterterms by 



(4.5) 



(4.6) 



5. Consistency condition 

It is very important to note that the variational formulas introduced in the previous 
section are assumed to be valid for finite parameters m,gE, i.e., they are assumed to be 
valid beyond perturbation theory. Though it is difficult to derive the variational formulas 
from first principles, it is easy to check their consistency. Especially we will check their 
consistency with the RG eqs. 

We first examine the transformation of the left-hand sides of the variational formulas 
under the RG. From eqs. (Q, (^), and O we find 

^■(O,,(r0..A„(r„))^,, 



Y{X^,,+V,,,,{m,gE))[-Q-{0,Ari)...0, 
k=i ^ 



gE 



(5.1) 



{l + P^{gE))[- — {0,,{r,)...0,^{r^))m,g, 

m,gE 

fc=l 
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and 

d f d 



dl V dgE 



m,gE 



n , 
fc=l ^ 



m,gE 

n 

m,gE- 

Next we consider the transformation of the right-hand sides of ( [1.3|) . Prom eqs. 
and (^]^) we obtain 

d 
Jl 



^{0*^0.,iri)...OM)m,gE 



and 
d_ 
Jl 



= +r,,,,(m,r7i.)) {0*^0,Ari)...0,irk)...0,„irr,)U,gE 

k=i 

n 

- 2tv^ Hrm^^j{m, Qe) (Oil {ri)...Oj{rk)...0^^{rr,))m,gE 

k=l 

^ f d \ 

+ XI ( Jf"^^^^^ + (1 + M9E))Crmk,j + [Cm,X + T{m, gE)],^j j 
k=l ^ ^ 

■ {0,,{n)...Oj{rk)...0,^{rn))m,gE 



m,gE 



^Y.{X,^,j+r,^,j{m,gE)) {0*EO,M)-Oj{rk)...0,^{rr,))m,gE 

k=l 

- PE\9E){0l,0,,in)...0,^{rr,)U,gE - rnPrn\9E){0*^0,Ar^ 

n 

-27v^Y.^E^'''J^'^^9E){0,Ari)...Oj{rk)...OM)m,gE 

k=l 

^ f d \ 
+ XI ( Jf^^'"^ + ^E\gE)cEi^,3 + ml3m\gE)cmi^,j + [ce, X + r(m, gE)],^^j J 

■ {0,Ari)-0,{rk)...0,Sr^))m,gE- 
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Here we get the He terms from the boundary integrals at r = 1, and 2n'^ is the volume 
integral of a unit 3-sphere. 



By comparing eq. (|5.1j ) and eq. (|5.3| ), and eq. (|5.2|) and eq. ( pT^) we obtain 

27i'^Hmim,gE) = d^T{m,gE) + ^Cmim.gE) 

dl 

+ (1 + f3mi9E))cm{'m, gE) + [c^(m, gE)-,X + r(m, gE)] 
27r^i7£;(m, gE) = dgj,V{m, gE) + ^^CE{m, gE) ^^'^^ 

+ PE'{9E)cEim, gE) + m(3rn'{9E)cmim, gE) 

+ [c£;(m, gE), X + r(m, gE)]. 
These formulas give relations between the operator coefficients and anomalous dimensions. 
We can rewrite these in a more suggestive way as follows: 

= + [Cm, $] + (^E{dgECm " dmCE + [ce, Cm]) 

(5.6) 

2'k'^He = dgj,^ + [ce,^] -m{l+ /3m){dgj,Cm - dmCE + [cE.Cm]), 

where 

$(m, gE)=X + V + (3eCe + m(l + /3^)c^. (5.7) 

Eqs. ( |5.6|) constitute a main result of this paper, and they are analogous to the relation 
between the anomalous dimensions and 1/e poles in the dimensional regularization with 
the minimal subtraction. 

We will come back to the geometric meaning of the expressions ( p.6| ) in sect. 7. Note 
that the anomalous 

dimensions F contribute only to the maximal parts of the coefficient functions. 
Now, by substituting eqs. ( p75|) into eqs. (|3.9| ) we can express the coefficient functions 
for an arbitrary r as total derivatives: 

Cm{r;m,gE) 

1 d 



2n'^mr'^ dlnr 
CE{r;m,gE) 

1 d 



G{r;m,gE){X + m(lnr)c^(m(lnr), 5f£;(lnr)) IG ^{r;m,gE 



G{r;m,gE){T{m{lnr),gE{lnr)) - Prn{gE)X 



2tt'^ 13 E{gE)r^ 91nr 

+ 13 E [qe (In r))cE (rn(ln r) , (je (In r) ) 

+ m(lnr)(/3^(^£;(lnr)) - /3^(5f^;))c^(m(lnr), ^^;(lnr)) ^^"^(r; m, ^i^;)^ . 

(5.8) 

11 



To derive this result we need to use eqs. ( p.l2|) , (|3.8|), and 



[X, r(m, qe)] = mdmTim, qe), (5.9) 

where the last equation is a consequence of eq. ( p.7|) . 

We can now integrate the above coefficient functions to obtain an expression for the 
subtractions Am and Ae, defined by ( [4.4|) , in terms of the anomalous dimensions F and 
the counterterms Cju^ce- We find 

Am{e;m,gE) 

= — f - X + ^(e; m, gE)XG~^{e; m, qe) 
m V 

+ G{t] m, 5f£;)m(lne)c^(m(lne), gE{\nt))G~^{e] m, qe)^ , 
AE{e;m,gE) 

if, , (5-10) 

- GXG-') - r(m, + G'r(m(lne),^E(lne))G'-i 



+ ( (In e) )c£; (m(ln e) , (In e) ) 



+ (/3m(fi's(lne)) -/3,„(5f£;))m(lne)c^(m(lne),5f£;(lne))jG' ^ j . 

Using the above formulas we can count the powers of 1/e in the subtractions. First 
we note, from (|3.6|) and (|3.7|) , that 



iX - GXG-^)i = 0, if Xi = 

(5.11) 

oc m"^^ , if Xi > Xj. 

Since G{l;m,gE) = 1, the term {X — GXG~^)ij corresponds to purely logarithmic con- 
tributions in e that vanish at Ine = 0. Similarly, we find that 

(T{m,gE) ~GT{m{\ne),gE{lne))G-^) oc m^^"^^' (5.12) 

are purely logarithmic in e without powers of e. Secondly we note, from ( [4.5|) , that 

(G'm(lne)c^(m(lne),^^(lne))G-i) ~ ,v(^i,(lne)) 

{GcEimilne),gEilne))G-')^^ ~ cg^^.,(^s(lne))^ 

12 



where Xi > Xi' > Xj/ > Xj. Especially we find 

(Gm{\ne)cm{m{\ne),gE{lne))G-A ~ c^i/,jv(^^(lne))m''^-'^^' 

V / i.i 



(5.14) 



These are not necessarily zero at e = 1 and involve terms 

that are finite as e ^ 0. In conclusion the anomalous dimensions F give only log- 
arithmic contributions to the subtractions Am,AE, while the non-maximal parts of the 
counterterms Cm.CE give contributions that are proportional to integral powers of 1/e, 
and the maximal parts Cm,CE give both logarithmic and finite contributions in e. In the 
absence of maximal counterterms there is no 

finite subtraction. Thus, the necessary subtractions are minimal when the maximal 
parts of the counterterms vanish. 



6. Commutativity 



In the previous section we checked the consistency of the variational formulas ( [4.3| ) 
with the RG equations. In this section we will perform another consistency check. The 
variational formulas give single derivatives of correlation functions with respect to m,gE, 

and we can use the formulas recursively to realize multiple derivatives in terms of 
multiple integrals over operators. 

We define two insertions {0'^0^...)m,gE (in this particular order) such that it equals 
—dgj^ {0^...)m,gE'i namely we apply the definition of to the correlation function in the 
definition of {0^...)m,gE- 



r — ri I >e 



{0*EO:^0,Ari)...0,^irn))n.,gE^lim 

e—^0 

n 

+ I]^m^„,(e;m,^7^;)(0^0,,(rl)...0,(rfc)...0,„(r,))^, 

k=i 

n 

-^dg^Arm^j{e; m, £?e)(C,, (ri)...Cj (rfc)...C,„ (r„))^, 



d\ {0*EOm{r)0,M)■■■O^Arn)) 



c 

m,gE 



k=l 



(6.1) 
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{Ol,OE{T')0,ATr)...0,Srn))'in,,, 

r' —ri\>e 



Similarly, we define {0^0*^0...)m,gE such that it equals —dm{0%0...)m,gE 
{0*^0%Oi,{r^)...OM)m,gE = Jim 

n 

k=l 

n 

-^dmAEii„j{e;m,gE){Oi^{ri)...Oj{rk)...Oi^ {rn))m, 



,9e 

k=l 



(6.2) 

In a similar way we can define {0^0^...)rn,gE and (C|;0|;...)m,gB- 

The insertions and O*^ are operator realizations of the partial derivatives —dm 
and —dg^ . Since the partial derivatives commute with one another, we must satisfy 

{0*EO*mOM...OM)m,gE = {0*mO%OM...OM)m,gE- (6-3) 

The evaluation of the difference O'^O^ — O^O'^ is straightforward but tedious, 

and we will sketch the calculation in Appendix A. For simplicity we give the result 
only for the case n — 1 (For the case of a general n, see Appendix A.): 



(6.4) 



{{0%0*^-0*^0%)Omm,gE 

— {^Em)i,j i^j) m,gE ~ {^gE'^rn dmC-E ~l~ [cs? C^] )i,j (Cj)m,gB • 

Here we define the curvature 

{^Em)i,j{^3)m,gE 

= / dVF.P. / d^r' {Om{r){OE{r')Om -CEiAr'-,m,gE)Oj{Q)) (6.5) 

J\>r Jl>r' 

- OE{r){Om{r')OiiO) - CmiAr';m,gE)Oj{0))) 

m,gEi 

where F.P. denotes an integrable part with respect to r. 
Thus, the consistency condition becomes 

dgECm - dmCE + [Cfi, Cm] = ^Em- (6.6) 

The left-hand side is reminiscent of a gauge field strength with Cm,CE as the gauge fields. 
We will elaborate on this view-point in the next section. It is important to emphasize that 
the curvature 0,Em is related essentially to a product of three operators close together, 
and it cannot be deduced from the OPE coefficients which have to do with the products 
of only two operators close together. 
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Before we dose this section we must ask if 

we get further consistency conditions from more 

number of insertions of and C^. Let us consider insertion of K Oj^'s and L O^'s. 
These are defined using ( [4 .31) recursively. We can classify all possible insertions into two 
classes, those of type 1 

and those of type 2 

where O* ... corresponds to the remaining K + L — 1 insertions. Suppose that the insertion 
of K' Oj^'s and L' O'^'s is independent of ordering if K' + L' < K + L. Now, using the 
variational formulas ( [4.3| ) only once, we get 



{0:^0*...0,Ari)...O.Arn))n.,,, = -d^{0\..0,Ari)...OM)m,o. 

Hence, by the assumption all correlation functions are equal in each class. The equality 
between the two classes amounts to show 

{0*^0],0\..0,Ar,)...0,^{r^))^,g, = {0],Ol,0\..0,,{n)...OM)^,g, (6.8) 

where 0*... stands for the remaining i^'+L— 2 insertions. Since {O* ...Oi^{ri)...Oi^{rn))m,gE 
is given as an integral over iiT + L — 2 + n-point correlation function (ignoring counterterms 
for simplicity), we can prove the equality by going through the same proof as given in Ap- 
pendix A for the case of two insertions. Therefore, we will not get any further consistency 
condition from multiple insertions. 



7. Unique convention for composite operators 

We would like to redefine the composite operators in such a way that the corresponding 
counterterms have vanishing maximal parts. The reason is that this convention is the 
simplest in which finite subtractions are unnecessary. 

The most general redefinition of the operators {Oi} that we can think of is given by 

0, = Y,N^,M.gE)0,, (7.1) 
j 
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where the matrix is invertible. The element Ni^j can 

be nonvanishing only if Xi~Xj is a non-negative integer, and Nij must be proportional 

to rn^^~^^ : 

dmN=-[X,N]. (7.2) 
m 

This condition is required in order to preserve the scaling dimensions of the redefined 
operators additively under the RG. With respect to qe, the element Nij must be a power 
series as required by locality of the theory. 

For the redefined operators O, the counterterms for 

the variational formulas (|4.3[) are given by 



= Nc^N-^ - d^N ■ 

(7.3) 

CE = NcEN-'-dg^N-N-\ 

Hence, the finite counterterms transform as gauge fields in the space of parameters m, qe- 
The inhomogeneous terms are maximal, and for the maximal parts alone we find 



(7.4) 

SB" 



Ie = NceN-^ -d,,N -N-^ 



The transformation properties ( [7.3|) imply that we can interpret [cj^^ce] as a gauge 
field on the two-dimensional theory space whose coordinates are m and qe- The operators 
{Oi} form a basis of an infinite dimensional vector bundle. The commutativity condition 

m 

gives the curvature of the gauge field in terms of an integral over a correlation function 
( |6.5| ) . Since the integral has no reason to vanish, the curvature is nonvanishing in general. 
The field $ defined by (|3) can be 

seen to transform as 

|. = iV$iV-i. (7.5) 

Thus, $ is a scalar field in the adjoint representation on the theory space. Eqs. ( |5.6| ) 
and ( |6.6|) give the vector field {Hm,HE) in terms of a covariant derivative of $ and the 
curvature fiEm- 

Hence, {Hm, He) is in the adjoint representation as expected. So much for geometric 
interpretations of our results. 
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The simplest convention for the finite counterterms is the vanishing maximal counter- 
terms: 

Cm = 0, Ce = 0. (7.6) 

But in general this is impossible, since the curvature ^Em, defined by ( |6.5|) , is nonvanishing. 
We remove the ambiguity of the convention by imposing 

CE{m,gE) = (7.7) 

and 

5m{m, QE = 0) = 0. (7.8) 

This is the analogue of the temporal gauge used in non-abelian gauge theories. For com- 
pleteness we will derive the matrix of transformation in Appendix B. The gauge con- 
ditions ( |7.7D and ( |7.8D fix up to a left multiplication by a constant matrix. Hence, the 
choice of the composite operators is uniquely determined up to a 
left multiplication by a constant matrix. 

Finally we ask whether the unique convention we obtained depends on the choice of 
the renormalized parameters m,gE- The most general redefinition of the parameters are 
given by 



m = mf{gE) 
9~E = gigE), 

where 

f{9E) = 1 + fi9E + y9e + - 

9{9e) =9e + -^9e + ■^9e + - • 
The corresponding changes in the partial derivatives are 

f) \ f) 

f{9E) 

d - = ^ d - m ^'^^^^ d 
9'{9e) 9'{9E)f{9E) 

The operators Om and Oe transform in the same way as a vector: 



(7.9) 



(7.10) 



(7.11) 



f{9E) 

Oe = — r Oe - m— — — - — - Om- 

9\9e) 9{9E)f{9E) 
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(7.12) 



This induces homogeneous transformations of the operator coefficients: 

1 



Cm{r;m,gE) 



f{9E) 



(7.13) 



g'ioE) 9'{9E)f{gE) 

From eqs. (|7.11|) and (|7.13|) we find that the counterterms also transform homogeneously 
as a vector: 

1 



Cm{m,gE) 
CE{rh,g~E) 



Cm{m,gE) 



f{9E) 

CEim.gE) - m ^ Cm{m,gE)- 



(7.14) 



9'{9e) 9'{9E)f{9E) 
Therefore, under the change of parameters the gauge condition ( [7.71) is not preserved, 



while the condition (|7.8|) is preserved. Thus, in order to remove the ambiguity in the choice 
of composite operators by the gauge conditions ( [f77|) and (|7.8|) , we need to introduce a 
convention for the choice of m and gE- This can be done by applying conditions analogous 
to (|7.7|) , (|7.8| ) to the conjugate operators Om-, Oe- 



8. Concluding remarks 

In this paper we obtained two main results. First we obtained the variational formulas 
( [1.3|) that realize the partial derivatives with respect to m, gE in terms of the operator 
insertions Oj^ , O*^ . Secondly we found a relation between the anomalous dimensions of 
composite operators and their operator product coefficients, given by eqs. (|5^). 

Using the second result we can compute the subtractions Am, Ae, which are necessary 
in the variational formulas (f4.3| ), in terms of the beta functions, anomalous dimensions F, 
and finite counterterms Cm.CE- 

The finite counterterms are not unique, since the composite operators are susceptible 
to redefinitions ( [7.1| ). Once we remove this ambiguity by imposing conditions ( |7.7| ) and 
( [f78|) , the finite counterterms become unique, and they are determined by eqs. ( ^?6|) and 

the commutativity condition (|6.6|). 

The non-maximal part of the finite counterterms corresponds to subtractions divergent 
like powers of the cutoff distance e. (See the last paragraph of sect. 5.) The dimensional 
regularization only sees logarithmic divergences and has no counterpart to the non- maximal 
finite counterterms. The maximal part of the finite counterterms gives finite subtractions 
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in the variational formulas. This results from our asymmetric treatment of conjugate 
operators when we evaluate higher order derivatives of correlation functions with respect 
to parameters. In the dimensional regularization each operator conjugate to a parameter 
is treated symmetrically. 

We would like to emphasize the conceptual aspects of our results rather than their 
potential usefulness for perturbative calculations. In fact some modifications will be neces- 
sary to apply the variational formulas ( |4.3| ) for perturbative calculations of the correlation 
functions of composite operators. The reason is that the operator Oe, conjugate to qe, 
roughly corresponds to 

~ ^ tr (F^,)2, (8.1) 
9e 

where the field strength is defined by 

F^. = d^A, - d^A^ + [A^, A,] (8.2) 
in terms of the gauge field, since the lagrangian is given by 

L = — tr (F^,)2. (8.3) 

9e 

Note that this is not an artifact of normalization. Even if we redefine a new A by ^JqeA, 
we still get Oe ~ F"^ /qe- The form of the operator Oe given by ( P7T| ) implies that the 
correlation functions involving in fact diverge in the limit qe 0. (But this does not 
upset our assumption of the analyticity of the operator coefficients C^, Ce with respect to 
qe- In eqs. ( p.l| ) we have taken the angular average, which eliminates these singularities.) 
This kind of complication is absent in the ((/>^)4 theory. Hence, we can apply the variational 
formulas (|4.3|) to evaluate multiple derivatives of the correlation functions of composite 

operators at A = 0. In this way we can relate the short distance singularities in the 
renormalized theory to the singularities we encounter in perturbative calculations. Then, 
we can interpret eqs. ( |5T6| ) (or ( |5.1C1| )) as a relation between the perturbative singularities 
and anomalous dimensions. Calculations of OPE coefficients for the conjugate operators 
have been done up to first order in A. 

I would like to thank Orlando Alvarez for discussions. 
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Appendix A. Calculation of 0*^0*^ - 0*^0*e 

In this appendix we will sketch the derivation of the commutator 



From the 



definition ( |0| ) we find, using (|4.3|) 



{d^dg^ - dg^d^){Oi{Qi))m, 



lim 



SB 



^ hm ( y|^,_|^^ {{Om{r')OE{r)-OE{r')0^{r))Om 



m,gE 



+ A^,,,(r7)(Oi,(r)0,(O))^,,,-As.,,(r7)(O^(r)O,(O)) 



— / (i^r Sg lim ( ditto ) 

Jl>r>e ^ ^ 

+ CE,,, lim ( / rfV (O^(r')O,(0))^,3^ + A^,- fc(e)(C>fc) 

- C^.,, hm ( / rfV {OE{r')Oj{0))m,gE + AEj,k{^){Ok) 



9E 



m.,gE 



m,gE ; • 



(A.l) 

Here Sg denotes the singular part with respect to the coordinate r. 

Due to the antisymmetry of the integrand under interchange of m and qe^ the first 
integral over r' can be restricted to r' < e. In the second integral over r', we must keep 
the entire range of integration, since the operation of taking the singular part spoils this 
antisymmetry. Hence, we obtain 



{drudge - dg^dra){Oi)ra, 



lim 



hm / rfV {OE{r){O^{r')O,{0)-CrmAr')Ojm 

r>e \ '^^^ ^ J£>r'>r] 



- OUr){OE{r')Om - CE^Ar')OM)) 



m,gE 



+ I d'r' {-CrmAr'){OE(r)0,m^,g, + C,,.,, (r') (0^(r)O,-(O)) 

'l>r'>e ^ 



m,gE 



d\ Sg lim I / d^r' {{OEir)Omir') - O^{r)OE{r'))O,{0))m,gE 

r>e '7^0 \Jr'>r, 

+ I d'r' (~CrmAr'){OEir)0,mm,gE +CE^Ar'){0Ur)0,m^,9. 

'l>r'>ri ^ 



+ {dgj^Cm - dmCE + [ce, Cm])i,j{Oj) 



m,gE ■ 



(A.2) 
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Now using operator product expansions we find 



lim 



d\ / {OE{r){Om{r')Om - CrmA^')0,m) 



m,gE 



(A.3) 



e>r' 



and 



Sg [ d^r' {0^{r')OE{r)Om)m,gE = I d^r' (O^(r')C^;,,,(r)O,(0))^,3, (A.4) 

Jr'>l Jr'>l 



and similar equations obtained by interchanging m and qe- 
Substituting the above into ( |A.2|) , we obtain 



hm 



d^r / rfV {OE{r){Om{r')Om - C^,,,(r')O,(0)) 

l>r><E ^Jl>r' 



- Om{r){OE{r')O,{0) - CB,,,(r')O,(0)))^, 



9E 



(A.5) 



- / rfV Sg (ditto) 

Jl>r>e ^ ' 



This gives O- 

In the above result the range of the double integral is restricted to r < 1, since we used 
r = 1 as the renormalization point when we introduced divergent subtractions Ae-,A^ in 
( [4.4|) . If we use a different renormalization point r = tq < 1 to define the subtractions 
Ae, Am, then we must change the finite counterterms CE,Cm to 



Ce = ce - I d r CE{r) 

'l>r>ro 



d^rCm{r). 



(A.6) 



l>r>ro 

Then eq. (|A.5| ) can be rewritten as 
i^r F.P. f d^r' 

•J ro>r J ro>r' 

+ [dgj^Cm — dmCE + [cE^'^m])^,j{^j)m,gE^ 



d'^r F.P. / rfV {OE{r){Om{r')Om - Cm^,,(r')O,(0)) 



(A.7) 
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where F.P. denotes a finite part with respect to r. We can take the radius tq as smaU as 
we want. This guarantees the locality of the double integral. 

In the presence of other composite operators the consistency condition (|A.5| ) is mod- 
ified to 



{d^dg^ - dg^dm){0,,{ri)...0^^{rn))m, 



SB 



k=i 



d^r F.P. / 

ro>|r— r-fel ^ ro>\r' -rk\ 



m,gE 



(A.8) 



where we can take vq smaller than the smallest distance between any two of the n points 



n, ...,rn. 



Appendix B. Explicit calculation of N{m^ qe) 

We have found that the counterterms c^jCe must satisfy the consistency condition 



dg^Cm -dmCE + [cE: Cm] = ^Er 



fB.ll 



where the curvature matrix flEm is defined by eq. ( ]6.5D . By taking the maximal part we 
obtain 



dgj^Cm — dfnCE + [cEj C^] — ^Er 



where the bars denote the maximal parts. 
Under a redefinition of the operators 



0[^N,Am,gE)0 



J' 
where 

Nij{m,gE) oc m'^^"'^^ 
the finite counterterms transform as 

c'm = NcmN-^ - dmN ■ N-^ 



c'e = NceN-^ - dg^N ■ N-\ 



(B.2) 



(B.3) 
(B.4) 

(B.5) 
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The non-maximal terms transform homogeneously, while the maximal terms transform 
inhomogeneously : 



dE = NcEN-'-dg^N-N-\ 
We wish to choose the matrix N{m, qe) such that 

dE{m,gE) = 

dmim, 0) = 0. 
The gauge condition ( p.7| ) is solved by 



(B.6) 



(B.7) 



(B.8) 



N{m, qe) = M{m) ■ TZexp 



9E 



dx CE(m, x) 



(B.9) 



where TZ implies the increasing ordering of x toward right, and M{m) is an arbitrary 
invertible matrix dependent only on m. The exponential is analytic in both m and qe- 
Since 

7V(m,0) = M(m), (B.IO) 



we find 

dmim, 0) = Mc,^(m, 0)M-^ - dmM ■ M'^. 
Hence, to satisfy the other gauge condition ( [B.8| ), we must choose 



(B.ll) 



M{m)=L- 7^exp 



dy Cm{y,Q) 



where L is an arbitrary invertible constant matrix. Note that M{m) is analytic in m. 
Finally we note that under the gauge conditions ( p.7|) , (|B.8|) , we get a relation 



c'm{^^9E)= dxn'Em{m,x) 

Jo 



[B.13) 
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